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Planck evidence for a closed Universe anda 
possible crisis for cosmology 


Eleonora Di Valentino’, Alessandro Melchiorri®2* and Joseph Silk?*° 


The recent Planck Legacy 2018 release has confirmed the presence of an enhanced lensing amplitude in cosmic microwave 
background power spectra compared with that predicted in the standard A cold dark matter model, where A is the cosmological 
constant. A closed Universe can provide a physical explanation for this effect, with the Planck cosmic microwave background 
spectra now preferring a positive curvature at more than the 99% confidence level. Here, we further investigate the evidence 
for a closed Universe from Planck, showing that positive curvature naturally explains the anomalous lensing amplitude, and 
demonstrating that it also removes a well-known tension in the Planck dataset concerning the values of cosmological param- 
eters derived at different angular scales. We show that since the Planck power spectra prefer a closed Universe, discordances 
higher than generally estimated arise for most of the local cosmological observables, including baryon acoustic oscillations. The 
assumption of a flat Universe could therefore mask a cosmological crisis where disparate observed properties of the Universe 
appear to be mutually inconsistent. Future measurements are needed to clarify whether the observed discordances are due to 


undetected systematics, or to new physics or simply are a statistical fluctuation. 


of cosmic microwave background (CMB) anisotropies has 

reported some unexpected results, revealing the possibility 
for new physics beyond the standard A cold dark matter (ACDM) 
model, where A is the cosmological constant'”. Indeed, while the 
inflationary predictions for coherent acoustic oscillations have been 
fully confirmed, a preference for a higher lensing amplitude, A,.,,; 
than predicted in the base ACDM model at about 3 standard devia- 
tions has been found in the temperature and polarization angular 
spectra. We argue that the An anomaly has profound implica- 
tions for some extensions to the ACDM model, such as the curva- 
ture of the Universe. The constraints from the PL18 CMB spectra 
on curvature, parameterized through the energy density param- 
eter Q,, are indeed surprising, suggesting a closed Universe at 3.4 
standard deviations (—0.007 > 2,>—0.095 at the 99% confidence 
level (CL)'*). 

As is well known, inflation theory naturally predicts a flat 
Universe*’. However, inflationary models with Q,<0 (refs. °*) are 
relatively simple to build, with primordial homogeneity and isot- 
ropy easier to achieve than in open models. An issue for closed 
inflation models is that to obtain Qg% —0.1, fine-tuning at a level 
of a few per cent is needed’. This does not sound very compelling, 
but it may still be acceptable, given the presence of a far more finely 
tuned cosmological constant. Closed models could also lead to a 
large-scale cut-off in the primordial density fluctuations, around 
the curvature scale R.=(c/H,)|Q,|°>1Gpce (where R, is the cur- 
vature radius, c is the velocity of light and H, is the current value 
of the Hubble constant), in agreement with the observed low CMB 
anisotropy quadrupole”. Confirmation of a positive spatial curva- 
ture would also have several implications for inflationary theory 
and, for example, severely challenge models of eternal inflation'®"'. 

Here, we show that, if indeed credible, the Planck preference for 
a closed Universe introduces a new problem for modern cosmology. 


| he recent Planck Legacy 2018 (PL18) release of observations 


Indeed, many of the current tight constraints on cosmological 
parameters are obtained by combining complementary datasets. A 
basic assumption in this procedure is that these different datasets 
must be consistent, that is, they must plausibly arise from the same 
cosmological model. Currently, two major experimental datasets are 
in tension with Planck: the determination of the Hubble constant 
by Riess et al.” is discrepant at the level of ~3 standard deviations 
(but see also ref. °), and the observations of cosmic shear by the 
Kilo Degree Survey 450 (KiDS-450) disagree at ~2 standard devia- 
tions'*’. Furthermore, the value of An, derived from the Planck 
lensing-generated four-point correlation function is consistent with 
the expectations of the ACDM model and in tension with the PL18 
power spectra™s. 

Although most of the remaining cosmological observables are 
considered to be in good agreement with PL18, these inconsisten- 
cies have already motivated several studies that have attempted to 
critically reassess the level of discordance’ or to resolve it with the 
introduction of new physics". 

The level of accordance between cosmological observables has 
hitherto been thoroughly investigated under the assumption of a 
flat Universe. We show here that when curvature is allowed to vary 
(as suggested by the PL18 CMB spectra), the statistical significance 
of the known tensions with PL18 increases and other discrepancies 
arise with several ‘local’ (at redshift z< 3) observables. The assump- 
tion of a flat Universe could, therefore, mask a cosmological crisis 
where disparate observed properties of the Universe appear to be 
mutually inconsistent. 

Before evaluating the tensions of the PL18 results with indepen- 
dent cosmological observables, we first check whether the PL18 
power spectra can provide an unbiased and reliable estimate of the 
curvature of the Universe. This may not be the case, since ‘geomet- 
rical degeneracy’ is present between cosmological parameters*”” 
For example, assuming the same inflationary parameters and 
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Fig. 1 | Preference for a closed Universe, 2, < 0, from Planck. Posterior 
(P) distributions normalized at the maximum value (Pra) for Q, from PL18 
temperature and polarization-simulated angular power spectra (assuming 
a fiducial flat ACDM model) and PL18 real data, adopting the baseline ‘Plik’ 
Planck likelihood and the alternative ‘CamSpec’ likelihood, respectively. 
For comparison, the posterior from the previous PL15** data release 

is also shown. 


reionization process, a flat cosmological model with matter den- 
sity 2,,=0.35, cosmological constant density @,=0.65 and 
H,=65kms"! Mpc"! produces an identical structure of the CMB 
angular spectrum at subdegree angular scales for a closed model with 
Q,,=1, Q,=0.15 (that is, Q,=—0.15) and H,=38.4km s~ Mpc™. 
Because of the form of the degeneracy, different closed models 
have identical CMB power spectra to that of a single flat model. 
The main consequence is that, after marginalization over the nui- 
sance parameters, the posterior on Qx is generally skewed towards 
closed models”. 

The situation changes with precise CMB measurements at arc 
minute angular scales: here, indeed, additional anisotropies induced 
by gravitational lensing are not negligible. Since gravitational lens- 
ing depends on the matter density, its detection breaks the geomet- 
rical degeneracy. The Planck experiment, with its improved angular 
resolution, therefore offers the opportunity of a precise measure- 
ment of curvature from a single CMB experiment. 

To confirm this hypothesis, we generated a Monte Carlo Markov 
chain (MCMC) analysis over simulated Planck (temperature and 
polarization) data, assuming the best-fit flat ACDM model and 
experimental noise properties similar to those presented in the 
PL18 release’. As we can see from Fig. 1, the expected posterior is 
centred around 2,=0, with a bound of Q,=0.00 +0.02 at the 68% 
CL. Potentially, an experiment such as Planck could constrain cur- 
vature with ~2% uncertainty, without any substantial bias towards 
closed models. 

For comparison, in Fig. 1, we have plotted the posterior from the 
PL18 real temperature and polarization power spectra, assuming 
the baseline Planck likelihood (see ref. *). As we can see, the pos- 
terior is reasonably centred on a closed model around Q,=—0.04. 
Integrating this posterior distribution over Qx we find that Planck 
favours a closed Universe (Q,<0) with 99.985% probability. 
Moreover, a closed Universe with Q= —0.0438 provides a bet- 
ter fit to PL18 with respect to a flat model, with a x? difference of 
Ayr ~ —11, where 72, is the effective best-fit chi squared from the 
MCMC chains’. 

This qualitatively shows the PL18 preference for a closed 
Universe, but does not statistically weight the additional param- 
eter (Qx). To better quantify the preference for a closed model, we 
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Fig. 2 | Degeneracy between curvature and lensing. Constraints at 68% 
(dark blue) and 95% (light blue) CLs in the Aiens versus 2x plane from PL18 
temperature and polarization data. A degeneracy between curvature and 
the Ans parameter is clearly present. Note that a model with <0 is 
slightly preferred with respect to a flat model with A,.,,,>1. 


adopt the deviance information criterion (DIC)”*', which takes 
into account the Bayesian complexity, that is, the effective number 
of parameters, of the extended model” and is defined as 


DIC =27. =x (1) 


where the bar denotes a mean over the posterior distribution. This 
quantity can be easily computed. We restrict the analysis to models 
with curvature in the range —0.2 <Q, <0; that is, we neglect open 
models because they are both disfavoured from observations and 
more difficult to realize in an inflationary scenario. We find that 
the Planck data yield ADIC = —7.4; that is, a closed Universe with 
Q,=—0.0438 is preferred, with a probability ratio of about 1/41, 
with respect to a flat model. 

We also compute the Bayesian evidence ratio by making use of 
the Savage—Dickey density ratio*’*’*’. Assuming the Savage-Dickey 
density ratio, the Bayes factor B, can be written as 


p(Ox|d, Mı) 
5 See an a a 2 
i m(Ox|M)) Ox=0 2) 


where M, denotes the model with curvature, p(Q,|d, M,) is the pos- 
terior for Q, in this theoretical framework, computed from a spe- 
cific dataset d, and 2(Q,|M,) is the prior on Q, that we assume to be 
flat in the range —0.2 < Qg <0. 

Applying the Savage-Dickey method to the Planck temperature 
and polarization, we obtain the Bayes ratio of 


| In Bo: | = 3.3 (3) 


that is, by assuming the so-called Jeffreys’s scale, we obtain strong 
evidence for closed models with Qx in the prior range [—0.2,0]. 
While the assumption of a larger prior would lead to weaker evi- 
dence, the preference from the data for a closed Universe is clear. 
This evidence could come from an unidentified systematic in the 
Planck data. However, as we can also see from the posteriors in Fig. 1, 
the preference for a closed Universe increases as we move from 
the Planck 2015 (PL15)* data release to the current PL18 release. 
Moreover, even assuming a significantly different procedure for the 
likelihood analysis’, and using the alternative CamSpec approach 
instead of the baseline Planck likelihood, the preference for 
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Fig. 3 | Curvature and shift in parameters. Constraints at 68% CL on cosmological parameters derived from two different multipole ranges 

(2 < £ < 800 and 800<é < 2,500) of the PL18 temperature and polarization data, assuming either aACDM model (left) or a closed model (right). 
Polarization data at low multipoles (2 < £ < 30) are included in both cases. The difference in the parameter constraints, present in flat ACDM, 
disappears when assuming a model with Q,=—0.045. The parameters for each panel are non-dimensional with the exception of Hy which has units of 


4 4 : ; a 
km s7 Mpc". laas maximum l; fi, minimum /. 


curvature is reduced, but is still well above 2 standard deviations 
with Ox = —0.03779;53; at the 95% CL’. In the case of CamSpec, we 
find Q,<0 with a 99.85% probability. Although the indication for 
a closed Universe is less significant with CamSpec, it is still present, 
showing that our result is not due to differences between analysis 
methods. 

The preference for closure in the Planck data is strongly con- 
nected to the higher lensing amplitude. This is evident from the 
parameter degeneracy between Ans and Qg, as shown in Fig. 2, 
where we report the two-dimensional (2D) constraints at 68% and 
95% CLs on Ans and Q, from the PL18 temperature and polariza- 
tion data’. The dark matter content can indeed be greater in a closed 
Universe, leading to a larger lensing signal, solving the A),,, anomaly 
and providing a robust physical explanation. As we can see, when 
a closed model is considered, A,,,,, is in agreement with the expec- 
tation of Ans =1. The amplitude of the lensing signal in Planck 
temperature and polarization data is precisely what is expected in 
a closed Universe. It is interesting to note that a ACDM +, analy- 
sis provides a marginally better fit to the ACDM + Aen, analysis, by 
Ay’ =—1.6, because closed models better fit the low-multipole data. 

As discussed in ref. ”, assuming a flat ACDM model, the values of 
the cosmological parameters obtained from the PL15 temperature 
angular spectrum in the multipole range 2 < £< 800 are ‘shifted’ with 
respect to those derived from the same Planck data relative to multi- 
poles in the range 800 < £<2,500. This tension is also present in the 
PL18 release’, and the inclusion of the A,n, parameter removes this 
difference. A key point of our paper is that the addition of curvature 
also solves this tension: in Fig. 3, we show that in a closed Universe 
with Q,=—0.045, the cosmological parameters derived in the two 
different multipole ranges, from PL18 temperature and polarization 
data, are now fully compatible. 
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However, if the PL18 power spectra suggest a closed Universe, 
the remaining cosmological observables are in strong disagree- 
ment with this. Let us now compare the Planck constraints with 
those coming from local observables, starting with baryon acous- 
tic oscillations. We first consider a combination of measurements 
given by the Six-degree Field Galaxy Survey (6dFGS)**, Sloan 
Digital Sky Survey-Main Galaxy Sample (SDSS-MGS)* and 
Baryon Oscillation Spectroscopic Survey Data Release 12 (BOSS 
DR12)”** (hereafter, we refer to this dataset simply as BAO, where 
BAO stands for baryon acoustic oscillations), as adopted by the 
Planck Collaboration’. The combination of this BAO dataset and 
PL18 power spectra produces a strong constraint on curvature, 
with Ox = 0.000818-0038 at the 95% CL', in excellent agreement 
with a flat Universe. Given the significant change in the conclu- 
sions from Planck alone, it is reasonable to investigate whether 
the BAO dataset is actually consistent with PL18. The level of con- 
cordance between the Planck and BAO data, even from a quali- 
tative point of view, is immediately clear from Fig. 4, where we 
plot the acoustic-scale distance ratio, Dy(z)/rj,. as a function of 
redshift, taken from several recent BAO surveys, and divided by 
the mean acoustic-scale ratio obtained by Planck temperature and 
polarization data adopting a ACDM + 2x model. Here, farag is the 
comoving size of the sound horizon at the time of the end of the 
baryon drag epoch, and D,, the dilation scale, is a combination 
of the Hubble parameter H(z) and the comoving angular diam- 
eter distance D(z): Dy(z) = (czDi,(z)/H (z))'/3. As we see, there 
is a striking disagreement between the PL18 power spectra and 
the BAO data. This can also be seen in Table 1, where we report 
the constraints on D,, and H(z) from the recent analysis of BOSS 
DR12 data” and the corresponding constraints obtained indirectly 
from Planck, assuming a ACDM model with curvature. Each of 
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Fig. 4 | Tension with BAO. Acoustic-scale distance (D,) measurements 
divided by the corresponding mean distance ratio from PL18 temperature 
and polarization power spectra in a ACDM + 2x model. The green bands 
show the 68% (dark green) and 95% (light green) confidence ranges. 
The data points correspond to the measurements at the 68% CL from the 
following experiments: 6dFGS*°, SDSS-MGS7” and BOSS DR12*° (the BAO 
dataset considered in this paper). We also report measurements from 
WiggleZ°’, the DES*°*, the DR14 Luminous Red Galaxy sample’, SDSS 
Quasars™ and BOSS Lya”. 


Table 1 | Comoving angular diameter distances D,, (corrected by 
the fiducial comoving sound horizon r, at the baryon drag epoch 
assumed in the survey's analysis (r4fa)) and Hubble parameter 
measurements from recent BAO observations from BOSS 
DR12*° compared with the corresponding quantities derived 
from PL18 power spectra assuming aA CDM + Q, model 


Observable Redshift BAO Planck Tension 
(68% CL) (68% CL) 

DuCastia/ ra) (Mpc) 0.38 1,518+22.8 1,843+100 2.96 
DuCastia/ ra) (Mpc) 0.51 1,977+26.9 2,3614115 3.00 
Dulrastia/ ra) (Mpc) 0.61 2,283+32.3 2,726+130 3.30 
Hlragia/ fa) 0.38 81.5+1.9 71.6+3.3 2.60 
(kms Mpc”) 

(ra sia/Ta) 0.51 90.5 +1.97 78.9+31 31o 
(kms Mpc™') 

Hlrasia/ fa) 0.61 97.3421 85.0+3.0 3.36 
(kms! Mpc’) 


the BOSS DR12 data points is in disagreement by about 3 standard 
deviations with the Planck power spectra. 

As we can see from Table 2, the PL18 72, best fit is worse by 
Ay’ = 16.9 when the BAO data are included? under the assumption 
of curvature. This is a significantly larger Ay’ than obtained for the 
case of ACDM (Ay’# 6.15). The BAO dataset that we adopted con- 
sists of two independent measurements (6dFGS* and SDSS-MGS”*) 
with relatively large error bars (that is, with low statistical weight; 
see Fig. 4), and six correlated measurements from BOSS DR12*. It 
is therefore not straightforward to determine the number of inde- 
pendent data points present in the BAO dataset and to estimate 
the disagreement between the datasets from a simple y’ analysis. 
Although several statistical methods have been proposed to quan- 
tify the discrepancy between two cosmological datasets", here, we 
check for consistency between two independent datasets, D, and D,, 
by evaluating the following quantity based on the DIC approach'*”: 


T(D,, D2) = exp{—F (D1, D2)/2} (4) 


Table 2 | Tensions between PL18 and BAO and CMB lensing 


Additional dataset Aven AN ata logioZ 
flat ACDM 
+BAO +615 8 0.2 
+ CMB lensing +8.9 9 0.6 
ACDM+Q, 
+BAO +16.9 8 -1.8 
+ CMB lensing +16.9 9 -0.84 


In the second column, we report the best-fit Ang with respect to the PL18 dataset alone. In the 
third column is the number AN,,,. of (correlated) experimental data points from the additional 
dataset. In the fourth column is the value of log o£ that quantifies the tension (substantial if less 
than —0.5, strong if less than —1.0). We note good agreement between the datasets in the case of 
flat ACDM. In contrast, statistically significant tensions arise when curvature is considered. 


where 


where DIC(D, UD,) is the DIC obtained from the combined analy- 
sis of the two datasets. 

Following the Jeffreys’s scale, the agreement/disagreement is 
considered ‘substantial if | log,, Z| >0.5, ‘strong’ if |log,, Z| >1.0 
and ‘decisive’ if | logio Z| >2.0. When log, Z is positive, then two 
datasets are in agreement, whereas they are in tension if this param- 
eter is negative. We show in Table 2 the values of log,, Z computed 
for the PL18 (D,) and BAO (D,) datasets in the case of ACDM and 
ACDM + Qx: For the ACDM model, there is reasonable agreement 
between the datasets (log,, Z = 0.2), but evaluating models with 
curvature results in strong disagreement log,, Z = —1.8) between 
Planck and BAO data. 

A second tension is present between PL18 power spectra and 
the constraints on the lensing potential derived from the four-point 
function of Planck CMB maps” (hereafter, called CMB lensing). 
Indeed, as discussed previously, the preference for Q,<0 in PL18 
is mostly due to the anomalous lensing amplitude at small angular 
scales’. However, this greater lensing amplitude is not seen in the 
CMB lensing data, which is consistent with a flat ACDM model. 
This can be seen in Fig. 5, where we compare the lensing-poten- 
tial power spectra best fits from the PL18 power spectra, obtained 
under the assumptions of curvature or flatness, with the CMB 
lensing data”. The flat ACDM model is in reasonable agreement 
with CMB lensing, while the PL18 best-fit Q,=—0.0438 model 
predicts a too large lensing amplitude (with the exception of two 
data points). As we can see from Fig. 4, the PL18 power spectra 
best-fit closed model predicts a lensing-potential spectrum that 
is very similar to the best fit obtained under ACDM+A,,,, with 
Ans = 1.191 (ref. °). 

A PL18+CMB lensing analysis yields Qx = 0.011793} at 
the 95% CL, bringing a flat Universe back into agreement within 
2 standard deviations, but still also suggesting preference for a 
closed Universe. However, it is interesting to quantify the discor- 
dance between PL18 and CMB lensing. As we can see in Table 2, the 
inclusion of CMB lensing in PL18 increases the best-fit chi squared 
by Ay’=16.9 in the case of ACDM+, (while in the case of the 
ACDM model, we have Ay’=8.9). The CMB lensing dataset con- 
sists of nine correlated data points. Even assuming these data points 
to be independent, the increase in y? when curvature is varied sug- 
gests there is tension at the 95% CL, while there is no significant 
tension in the case of flatness. Also in Table 2, we report the values 
ofthe Z quantity. We identify substantial agreement between PL18 
and CMB lensing in the case of a flat Universe (log,, Z = 0.6) that 
changes to substantial discordance (log,, Z = —0.84) when curva- 
ture is allowed to vary. 
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Fig. 5 | Tension with CMB lensing. The solid lines are the theoretical 


predictions for the best-fits CMB lensing power spectra (a) from the 
PL18 angular spectra in the case of ACDM, ACDM + Qx and ACDM + Aiens 
models, respectively. The ACDM+Q, model has 2,=—0.0438, while the 
flat ACDM + Ans model has A,,,,;=1.191. The grey bands (errors at 68% CL) 
are the CMB lensing conservative experimental band powers 

extracted from the PL18 trispectrum data, while the black dots are 

the measured values. 


In conclusion, if the assumption of a flat Universe is removed 
and curvature is permitted, as preferred by the PL18 power spectra, 
we find strong disagreement between Planck and BAO data, and 
substantial disagreement between Planck and CMB lensing data. 

It is interesting to investigate whether astrophysical measure- 
ments that are already in tension with Planck under the assumption 
of a flat Universe are still in disagreement when curvature is con- 
sidered. In a ACDM+, model, PL18 power spectra provide the 
constraint: Hp = 54.473} at the 68% CL’. This is now in tension at 
the level of 5.2 standard deviations with respect to the conservative 
constraint of H, = 73.52 + 1.62 at the 68% CL from the Riess et al.” 
analysis (R18). The inclusion of curvature, therefore, significantly 
increases (by ~48%) the tension between Planck and R18. 

A similar increase in the tension is present with cosmic shear 
data from KiDS-450. In Fig. 6, we show the 2D constraints in the 
Ox the amplitude of mass fluctuations, versus 2, plane from KiDS- 
450 and PL18 power spectra under the assumption of curvature. 
For comparison, we also include the Planck constraint under flat 
ACDM (the KiDS-450 bound is just slightly different when flat- 
ness is assumed). As we can see, there is a significant shift in 
this plane for the PL18 constraint when moving from a flat to a 
closed Universe, which increases the discrepancy with KiDS-450. 
Considering S; = 0,(@,,/0.3)°°, from the PL18 power spectra, we 
obtain S¿=0.981 + 0.049 at the 68% CL, a value that is now about 
3.8 standard deviations from the KiDS-450 result. Cosmic shear 
measurements have also recently been made by the Dark Energy 
Survey (DES)* and by the Subaru Hyper Suprime-Cam (HSC). 
These measurements are reasonably consistent with the PL18 result 
in a flat Universe. However, assuming that the reported constraint 
on the S, parameter depends weakly on Q,, we find that once curva- 
ture is allowed, the PL18-derived determination of S, is discordant 
at more than 3.5 standard deviations with the DES and at more than 
3 standard deviations with the HSC. In practice, when curvature is 
included, not only the significance of the tension with KiDS-450 
increases, but PL18 is now also significantly discordant with recent 
cosmic shear surveys such as the DES and HSC. 

Until now, we have studied the compatibility of single datasets 
with PL18. However, analyses are usually performed by combining 
multiple datasets. It is therefore interesting to address the compat- 
ibility of Planck with combined datasets. In Fig. 7, we show the con- 
fidence region at the 95% CL from a BAO +SNe-Ia+ BBN dataset 
(where SNe-Ia refers to type-Ia supernovae and BBN to big bang 
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Fig. 6 | Tension with cosmic shear measurements. Discordance between 
PL18 and the KiDS-450 cosmic shear survey in the c versus Q,, plane. 
Darker shaded areas show constraints at 68% CL while lighter shaded 
areas show constraints at 95% CL. The tension already present under 

flat ACDM (at about 2.3 standard deviations) is increased to more than 3.5 
standard deviations when curvature is incorporated into the analysis 

of PL18. 


nucleosynthesis) and the 68% and 95% CLs from the PL18 power 
spectra on the 2, versus H, plane. As we can see, there is strong 
tension between the Planck result and that from the combined 
BAO +SNe-Ia+BBN analysis. In principle, each dataset prefers 
a closed Universe, with the BAO+SNe-Ia+BBN dataset provid- 
ing just an upper limit of Q,<-—0.124 at the 68% CL. However, 
whereas the Planck result prefers Hy = 54733kms"!Mpc™ at 
the 68% CL, we find that the BAO +SNe-Ia+ BBN dataset gives 
Hy)=79.6+6.8kms7! Mpc"! at the 68% CL; that is, they are incon- 
sistent at the level of 3.4 standard deviations. Moreover, the 
BAO +SNe-Ia + BBN data prefer lower ages of the Universe (t,), with 
to = 11.7319 3? Gyr at the 68% CL, which is in modest tension with 
the recent age determinations (t.) of the stars 2MASS J18082002- 
5104378 B (ref. ®) and HD 140283, of t.= 13.535 + 0.002 Gyr and 
~=13.5+0.7 Gyr (refs. **“*), respectively. 

As we can see from Fig. 7, the probability contour plots from the 
BAO + SNe-Ia+ BBN analysis are rather broad. It is therefore inter- 
esting to include further observables to improve the constraints. We 
consider, separately, CMB lensing, the R18 determination of the 
Hubble constant and the observed angular size of the sound horizon 
at recombination, A, = 1.04116 + 0.00033 (where the subscript MC 
stands for Monte Carlo), in a ACDM + 2x model from PL18. We 
show the results of this kind of analysis in Fig. 8. The inclusion of 
the 0,;- prior from Planck shifts the constraints towards a flat ACDM 
model with Q,=0.0016+0.0075 at the 68% CL. The inclusion of the 
CMB lensing dataset also significantly improves the constraints, with 
Q,=0.00+0.01 at the 68% CL. Both BAO +SNe-Ia+ BBN + uc 
and BAO + SNe-Ia+ CMB lensing combinations provide evidence 
for a flat Universe, with good consistency between the datasets. We 
may argue that when deriving constraints under the assumption of a 
flat Universe, it would be more conservative to use these data combi- 
nations instead of PL18, since they are consistent with a flat ACDM 
model and do not show significant internal tensions. However, 
these data combinations still show a significant discordance with 
PL18 power spectra. Considering the parameter constraints derived 
from the BAO + SNe-Ia+ CMB lensing dataset, we indeed find dis- 
agreement with PL18 at 2.4 standard deviations in Q,, amounting 
to 2.7 standard deviations in H, and 2.9 standard deviations in S,. 
When we consider the combination of BAO + SNe-Ia+ BBN+ R18, 
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Fig. 7 | Tension with combined data. Contour plots at 68% (inner 

contour) and 95% (outer contour) CLs from PL18 and 95% CL region from 
BAO + SNe-la + BBN datasets in the Ho versus Q, plane and as a function of 
the age of the Universe. 


we find 2,=—0.091 + 0.037 at the 68% CL, which, again, provides 
an indication of a closed Universe (see Fig. 8). Both datasets pro- 
vide good best-fit chi-squared values, and it is impossible to dis- 
criminate one result over another from the statistical point of view. 
As we can see from Fig. 8, there is good agreement between the 
BAO + SNe-Ia + CMB lensing and BAO + SNe-Ia+ BBN + Oy, data- 
sets, while both are in significant tension at the level of 2.5 standard 
deviations with BAO + SNe-Ia + BBN +R18. 

In summary, the PL18 CMB power spectra provide a statisti- 
cally significant indication for a closed Universe. A closed Universe 
solves the internal tensions present in the Planck dataset on the 
value of the cosmological parameters derived at different angular 
scales. Positive curvature is also marginally suggested by the ages 
of the oldest stars (see, for example, refs. ****) and, in a combined 
analysis with the Aiens parameter, slightly favoured by the low CMB 
quadrupole. 

Apart from these arguments, none of the local cosmological 
observables currently favour a closed Universe, and most of them 
are consequently in significant discordance with PL18. BAO surveys 
disagree at more than 3 standard deviations. CMB lensing is in ten- 
sion at the 95% CL. The R18 constraint on the Hubble constant is in 
tension with PL18 at more than 5 standard deviations, while cosmic 
shear data disagree at more than 3 standard deviations. 

These inconsistencies between disparate observed properties of 
the Universe introduce a problem for modern cosmology: the flat 
ACDM model, de facto, does not seem any longer to provide a good 
candidate for concordance cosmology, given the PL18 power spec- 
tra preference for a closed model. At the same time, a closed model 
is strongly disfavoured by a large number of local observables. 

Clearly, a possible solution to this problem would be to speculate 
about the presence of hitherto undetected systematics in the PL18 
release. However, the statistical significance for a closed Universe 
increases when moving from PL15 to the PL18 release. We point 
out that the Wilkinson Microwave Anisotropy Probe satellite 
experiment”, after nine years of observations, also produced. the 
constraint Ox = —0.03770-*45 at the 68% CL, fully compatible with 
the Planck result. Finally, we have shown that discordance is also 
present between the R18 and the CMB lensing datasets once they 
are both combined with the BAO and SNe-Ia data. In practice, there 
is currently no supporting evidence that could lead us to believe 
that the observed inconsistencies are due to systematics in the PL18 
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Fig. 8 | Tensions in combined data. Contour plots at 68% and 

95% CLs (dark shading and light shading, respectively) from 

the BAO + SNe-la+ BBN + H,, BAO + SNe-la+ BBN + mc and 

BAO + SNe-la+ CMB lensing datasets in the Q, versus Ho plane. Also 
shown are the constraints on the Hubble constant from R18”. The black line 
is the mean value on Hubble constant from R18 and the yellow lines delimit 
the 68% CL region. 


data rather than in the low redshift measurements. Moreover, local 
probes are expected to be more contaminated by astrophysical sys- 
tematics and/or nonlinearities with respect to CMB anisotropies. 

If there are indeed no systematics in the Planck data, then the 
currently observed discordances may indicate the need for new 
physics and call for drastic changes in the ACDM scenario (see, for 
example, refs. 145-4), 

A third possible way is to consider the PL18 constraint on 2, 
as a, now reasonably unlikely, statistical fluctuation. Fortunately, 
future measurements will fully confirm or falsify current tensions 
and the PL18 evidence for curvature**’. In the meantime, we argue 
that the tensions with ACDM present in the PL18 release should 
not be discarded merely as a statistical fluctuation, but must be seri- 
ously investigated, since at face value, they point towards a drastic 
rethinking of the current cosmological concordance model. 


Methods 

Planck results. Most of the results presented in this paper have been obtained 
using the publicly accessible PL18 parameter chains, available at http://pla.esac.esa. 
int/pla/#cosmology. 


MCMC analysis. The results in Figs. 1-3 and 6-8 have been obtained by 
performing an MCMC analysis using the CosmoMC code”. All runs have reached 
a convergence such that R—1<0.02, where R is the Raftery-Lewis MCMC 
diagnostic function. For the Planck data, we considered the recently released PL18 
likelihood’. The 2D contour plot in Fig. 3 has been obtained assuming an eight- 
parameter model, where we assume the standard six parameters of the ACDM 
model (baryon density Q,h? (where h is the reduced Hubble’s constant), CDM 
density Q.h’, primordial amplitude A,, primordial spectral index n, optical depth 
to reionization t and sound horizon angular size 6,) plus Q, and Ans. The results 
in Figs. 7 and 8 are obtained from a BAO + SNe-Ia+ BBN dataset, where BAO is 
the compilation of baryon acoustic oscillations described in the paper, SNe-Ia are 
the luminosity distances of 1,048 type-Ia supernovae from the recent Pantheon 
catalogue”! and BBN is a prior on the baryon density of Q,h? = 0.0222 + 0.005 
derived from measurements of primordial deuterium” assuming big bang 
nucleosynthesis. 


Figure 3. The results in Fig. 3 have been obtained using the ‘clik change lrange 
bychannel pol.py’ tool, publicly available from the PL18 likelihood release’. 


Simulated Planck data. The simulated constraints on Qx in Fig. 1 have been 
obtained using a modified version of the makePerfectForecastDataset.py tool in 
CosmoMC to match the PL18 ACDM constraints. 
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Planck versus BAO comparison. Figure 4 has been obtained using a modified 
version of the BAO-data-z.py tool in CosmoMC. The log,,Z_ quantity has been 
computed either from the Planck parameter tables available at ref. ° or from a BAO- 
only MCMC run. 


Data availability 

The data that support the plots within this paper and other findings of this study 
are available at http://pla.esac.esa.int/pla/#cosmology or from the corresponding 
author upon reasonable request. 


Code availability 
All of the codes used to produce the presented results are publicly available. See 
Methods for more details. 
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